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ABSTRACT 


An  objective  analysis  technique  known  as  CRAM  (Conditional  Relaxation  Analysis 
Method)  has  been  programmed  for  the  IBM  7090  to  compute  ballistic  winds,  tempera¬ 
ture,  and  density  over  a  region  of  southeastern  Arizona  in  the  vicinity  of  Ft.  Huachuca. 
The  various  components  and  available  options  of  the  program  are  described.  A  set  of 
analyses  produced  from  initial  teste  is  shown  and  the  ballistic  quantities  derived  are 
compared  with  those  of  the  previous  quarterly  report. 

Corrections  based  on  a  single  station  (Ft.  Huachuca)  were  applied  to  concurrent 
artillery  firings  for  19  cases.  The  residual  errors  derived  from  these  cases  will  be 
used  as  a  control  for  later  experiments.  The  distribution  of  the  residual  errors 
relative  to  the  target  strongly  Suggests  a  bias  in  range  errors,  the  average  impact 
location  being  nearly  200  meters  beyond  the  target. 
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1.0  INTRODUCTION 


Work  in  the  third  quarter  centered  around  the  checkout  and  initial  test  runs  oi 
the  GRAM  (Conditional  Relaxation  Analysis  Method)  objective  analysis  program.  These 
tests  were  concerned  with  variations  of  forcing  functions,  initial  guesses,  relaxaU.n 
and  smoothing  procedures,  and  scaling  of  the  vertical  coordinate  as  applied  to  the  four 
analysis  parameters:  u-  and  v-wind  components,  temperature,  and  density. 

The  associated  gun  data  resulting  from  the  test  firings  of  a  pair  of  8-inch 
howitzers  are  being  studied.  A  small  computer  program  has  been  written  that  applies 
drift  and  coriolis  corrections  to  a  given  impact  point,  and  meteorological  corrections 
of  ballistic-wind  components,  temperature,  and  density.  The  program  was  run  on  a 
portion  of  the  January— February  1965  data  sample.  Results  are  described  in  Section 
3.0. 

The  third  activity  of  the  quarter  involved  the  preparation  of  a  post-processing 
program  that  will  operate  on  the  output  of  the  objective  analysis  program .  The  purpose 
of  this  new  program  is  to  compute  ballistic  corrections,  error  fields,  and  certain 
statistical  summaries  of  the  data. 


2.0  OBJECTIVE  ANALYSIS  EXPERIMENTS 
2.1  Description  of  GRAM  Program 

CRAM  is  the  principal  analysis  tool  for  this  study  to  provide  analyzed  maps  of 
winds,  temperature,  and  density.  The  main  components  of  the  CRAM  program  are: 

(a)  the  initial-guess  technique,  (b)  correction  procedures,  (c)  relaxation  methods, 

(d)  smoothing,  and  (e)  ballistic  computation.  The  description  of  these  components  that 
follows  points  out  the  features  of  this  IBM  7090  program. 

2.1.1  The  Initial  Guess 

CRAM  requires  the  use  of  an  initial-guess  field,  which  undergoes  subsequent 
corrections.  In  principle,  the  initial-guess  field  can  be  generated  in  a  variety  of  ways; 
Two  rather  reasonable  approaches  are  to  apply  a  surface-fitting  technique  to  the  data 
or  to  use  the  previous  analysis  (e.g.,  two  hours  earlier),  if  one  is  available,  or  some 
combination  or  “blend”  of  the  two  approaches.  Our  version  of  the  program  uses  these 
surface  fitter;  and  persistence  approaches. 

In  the.  surface  fitting  procedure,  the  three-dimensional  surface  is  sought  which 
best  fits  the  observed  field  of  data.  The  equation  is  of  the  form 

(p  =  +  a  x  +  a0x0  +  a0x„  +  ...  +  a  x  (2-1) 

0  11  22  33  pp 

where  the  dependent  variable  <p  is  the  variable  to  be  analyzed  over  an  evenly-spaced 
field  of.grid  points,  the  a’s  are  the  coefficients  and  the  x’s  are  the  independent  vari¬ 
ables  which  are  functions  of  the  locations  of  the  observations  of  <p,  i.e.,  x,  y,  and  z. 
Because  there  are  a  large  number  of  terms  that  can  be  generated  from  the  possible 
combinations  of  x,  y,  and  z,  many  of  which  may  not  be  significant  in  explaining  the 
variability  of  (p,  it  seems  reasonable  to  attempt  to  reduce  the  number  of  possible 
“predictors”  through  the  use  of  a  systematic,  stepwise  screening  procedure.  The 
method  chosen  is  that  of  screening  regression.  From  an  array  of  possible  predictors 
(x,  y,z,  xy,  xz,  etc.),  the  screening  procedure  first  selects  the  one  that  has  the  highest 
linear  correlation  with  the  prcdictand  in  question.  This  predictand  is  then  held  con¬ 
stant,  and  partial-correlation  coefficients  between  the  predictand  and  each  of  the 
remaining  predictors  are  examined;  the  predictor  now  associated  with  the  highest 
coefficient  is  the  second  one  selected.  Additional  predictors  are  chosen  similarly  until 
a  selected  predictor  fails  to  explain  a  significant  additional  percentage  of  the  remaining 
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variance  of  the  predicland.  The  coefficients^  (a^,  a^,  etc.)  are  then  derived  by  tlie 
method  of  least  squares.  The  initial-guess  field  is  then  generated  by  solving  the 
derived  equation  at  each  of  the  grid  points  in  the  analysis  volume. 

The  program  is  written  such  that  it  is  possible  to  compute  an  initial  guess  from 
a  blend  of  persistence  and  surface  fitting.  This  is  done  to  incorporate  continuity,  as 
the  observations  are  taken  at  two-hour  intervals  on  alternate  days  between  0600  and 
1400  MST.  For  the  first  observation  time  of  the  day  (0600  MST),  -however,  nonveight 
is  given  to  persistence  because  the  previous  available  observation  ismsually  forty 
hours  old.  The  formula  used  in  the' program  to  calculate  the  initial-guess  value  of  «p 
at  a  grid  point  is 


<P. 


I.G. 


W 

1  p 

- t U  +  - l-r-  <£ 

1+W  surface  fit,  1+W  persistence 


(2-2) 


where  W  =  weight  of  persistence.  By  assigning  a  value  of  zero  to  W^,  persistence 
is  given  no  weight  and  the  initial  guess  is  based  solely  on  the  surface  fit.  A  value  of 
one  for  W  assigns  half  the  weight  to  each  while  a  very  large  value  of  W,  essentially 
gives  all  the  weight  to  persistence. 


2.1.2  Correction  Procedures 

After  it  is  determined,  the  initial-guess  field  is  “corrected”.  For  each  observa¬ 
tion,  the  difference  between  the  observed  value  and  the  value  computed  for  that  loca¬ 
tion  by  interpolating  among  the  initial-guess  values  at  the  surrounding  grid  points  is 
computed.  This  difference  is  then  translated  to  the  nearest  grid  point  and  added  to 
that  grid  point  as  a  correction  to  the  initial  guess.  When  a  grid  point  is  subject  to 
multiple  corrections  (from  several  different  observations),  an  overall  correction  is 
computed  as  the  arithmetic  average  of  the  several  individual  corrections. 

2.1.3  Relaxation  Methods 

The  relaxation  process  is  the  next  step  in  the  objective  analysis  procedure.  Wb 
compute  for  all  grid  points,  except  those  on  the  boundaries  and  those  which  have  been 
corrected  by  data,  new  values  that  satisfy  the  Poisson  equation, 

v%(i,j,k)  =  F(i,j,k)  (2-3) 

where  <p  is  the  value  of  the  analysis  parameter  at  grid  point  (i,j,k),  F  is  a  forcing 
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function  which  may  define  some  specified  feature  of  the  field  of  <£,  and  “  is  the  finitc- 
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difference  Laplacian  operator  in  three  dimensions.  The  Laplacian  of  a  good  initial- 
guess  field  provides  a  good  forcing  function.  In  the  program,  the  formula  is  written 
as  .follows: 


^ <f>(i J,k)  =  $(i+l,j,k)  -i-  j,k)  +  $(i,j+l,k)  -  0(ij-l,k) 

£(ij,k)j 


-  4?(i,j,k}  4-  KA[0(ifj,k+l) 


-  K^M.k)  -  $(i,  j,k-l}]  (2-4) 

where  and  Kg  are  coefficients  pertaining  to  the  vertical  gradient  computed  above 
arid  below  the  analysis  grid  point,  respectively. 

The  horizontal  and  vertical  space  gradients  have  been  expressed  separately 
in  this  equation..  There  are  several  practical  reasons  for  doing  this,  all  having  to  do 
with  the  scaling  of  the  vertical  coordinate  used  in  the  analysis  technique.  The  basic 
horizontal  grid  length  is  10  km,  while  in  the  vertical  the  grid  length  is  about  an  order 
of  magnitude  smaller.  Measuring  gradients  over  this  smaller  distance  is  generally 
compensated  for  by  the  fact  that,  for  most  atmospheric  variables,  the  vertical  grad¬ 
ients  are  much  larger  than  the  horizontaLgradients  so  that  one  may  be  justified  in 
using  Ka  =  Kg  =  1  in  Eq.  (2-4).  The  relative  inequalities  between  the  grid  distances 
above  and  below  the  grid  point  which  arise  from  the  uneven  spacing  of  the  artillery 
zones  may  be  accounted  for  by  assigning  different  values  to  K^  and  Kg.  Also,  it  is 
possible  in  the  program  to  consider  the  objective  analysis  problem  in  a  two  dimen¬ 
sional  framework  by  simply  setting  K  =  K  =  0. 

A  13 

The  solution  of  Eq.  (2-3)  is  accomplished  by  relaxation  methods.  For  each  of 
the  grid  points'  where  new  values  are  being  computed,  the  residual,  R,  is  computed  by 

R(i,j,k)  =V20(iJ,k)  -  F(i,j,k).  (2-5) 

New  estimates  of  (p  are  then  computed  by 

?(i,j,k)  =Q-R(i,j,k)  +  <£(i,  j,  k)  (2-6) 

where  o,  a  relaxation  coefficient,  is  an  input  parameter.  New  residuals  are  then 
computed.  The  iterative  procedure  continues  until  all  residuals  arc  less  than  a 
specified  value  (e). 
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2.1.4  Smoothing 


A  smoothing  operator  has  been  developed  which  is  designed  to  eliminate  from  the 
analysis  undesirable  small-scale  features.  The  expression  used  is 

Sl$(i,j,k)]  =  W,  <?.{i,j,k)  *  \V0[$(i+l,j,k)  +  9(i-l,j,k) 


+  ^(i.  jrl.k)  +  ^(i,j-l,k)] 

+  W3[0(U,k-tt)  -  <?(i,j,k-l)]. 


(2-7) 


Bt 


B. 


W  =  - ■ - - - 

1  6  +  4B  -f2B 

\  h  v 


W  =  - - - \V  =  - : - 

2  6  -r  4B,  +2B  3  6  -s-  4B,  +2B 


v  / 


The  parameters  B^  and  B^  control  the  degree  of  smoothing  to  be  performed. 

A  value  of  zero  assigned  to  each  results  in  the  entire  weight  being  given  to  $(i,  j,k), 

i.e.,  no  smoothing.  Large  values  of  B^  and  B^  result  in  heavy  smoothing.  Also,  U  e 

smoothing  can  be  restricted  to.  two  dimensions  by  setting  B  =  0. 

v 

2.1.5  Ballistic  Computation 

The  CRAM  Program  produces  analyses  of  u-  and  v-wind  components,  tempera¬ 
ture,  and  density  for  each  of  ten  artillery  zones.  These  zones  are  listed  in  Table  2-1 
for  an  analysis  area  17  grid  units  in  the  I-direction  (from  1=50  to  1=66)  and  15  grid 
units  in  the  J-direction  (from  J=345  to  J=359)  or  about  170  km  x  150  km, 

A  portion  of  this  grid  is  shown  in  Fig.  2-1.  The  analysis  area  is  purposely  made 
large  so  that  analysis  problems  associated  with  the  boundaries* are  far  enough  removed 
to  not  contaminate  the  analysis  in  the  region  of  interest,  i.e.,  near  the  firing  area. 

The  ballistic  fields  are  produced  within  the  program  by  applying  the  zone  weighting 
factors  that  are  given  in  Table  2-1.  A  minor  modification  in  the  program  will  permit 
ballistic  computations  of  “lines”  other  than  merely  line  10. 

2.2  Analysis  Examples 

The  five  time  periods  of  25  January  1965  were  chosen  to  illustrate  the  CRAM 
Program.  This  is  the  same  data  sample  that  was  analyzed  by  hand  and  reported  in  the 
previous  quarterly  report*.  Results  of  these  analyses  should  be  considered  as  being 


*Ostby,  F.  P.,  1966:  Ballistic  Winds  Study,  Quarterly  Report  No.  2,  Contract  No. 
DA  28-043  AMC-01377(E),  The  Travelers  Research  Center,  Inc. 
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TABLE  2-1 

ZONE  STRUCTURE  AND  WEIGHTING  FACTORS  FOR  BALLISTIC 
COMPUTATIONS  (maximum  ordinate:  8,000  m) 


Zone 

Height  (m) 

Weighting  factors* 

Base 

Top 

Density 

1 

Surface 

200 

0.01 

0.00 

0.03 

2 

200 

500 

0.02 

0.01 

0.04 

3 

500 

1000 

0.02 

0.02 

0.07 

4 

1000 

1500 

0.04 

0.02 

0.07 

5 

1500 

2000 

0.03 

0.03 

0.07 

6 

2000 

3000 

0.07 

0.07 

0.13 

7 

3000 

4000 

0.08 

0.13 

0.12 

8 

4000 

5000 

0.09 

0.22 

0.11 

9 

5000 

6000 

0.09 

0.23 

0.11 

10 

6000 

8000 

0.55 

0.27 

0.25 

♦From  FM  6-16,  Tables  for  Artillery  Meteorology. 


of  a  preliminary  nature,  with  further  testing  of  the  program  options  being  required. 
The  analyses  are  shown  in  Figs.  2-2  through  2-11. 

The  initial-guess  procedure  of  blending  a  surface  fit  with  persistence,  described 
in  Section  2.1.1,  was  net  ready  in  time  for  this  particular  set  of  analyses.  An  alter¬ 
nate  procedure  with  a  four-dimensional  surface  fit  incorporating  time  (t)  was  used. 

The  screening  regression  technique  was  applied  to  parameters  generated,  from  various 
combinations  of  x,y,  z,  and  t.  The  resulting  set  of  equations  showed  considerable 
dependency  on  the  vertical  coordinate,  with  z  terms  being  selected  first  for  each  of  the 
four  equations. 

In  these  initial  tests  we  have  used  a  forcing  function  of  zero  [F  =  0  in  Eq.  (2-3)] 
and  have  also  set  =  0  in  Eq.  (2-4).  Testing  of  the  relaxation  coefficient,  a, 

indicates  that  the  number  of  iterations  are  minimized  with  a  =  0.22.  The  smoothing 
operator  was  employed  using  =  1.5  and  =  0,  which  give  half  the  weight  to  the 
grid  point  in  question  and  the  other  half  to  the  average  of  the  four  surrounding  grid 
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Fig.  2-1.  Objective  analysis  grid  showing  observing  stations.  Grid  units  are  based  on  the  10,  000-m 
Universal  Transverse  Mercator  Grid  (Zone  12),  with  the  last  four  digits  of  the  grid  numbers  omitted. 


Fig.  2-2.  Ballistic  wind  analysis  using  CRAM:  0600  MST,  25  January  1966 
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Fig.  2-3.  Ballistic  wind  analysis  using  CRAM:  0800  MST,  25  January  1965. 
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Fig.  2-4.  Ballistic  wind  analysis  using  CRAM:  1000  MST,  25  January  1965. 
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Fig.  2-5.  Ballistic  wind  analysis  using  CRAM:  1200  MST,  25  January  1965. 


Fig.  2-6.  Ballistic  wind  analysis  using  CRAM:  1400  MST,  25  January  1965. 


\-\i 


-r3  t.i  r/S  O  >  /jo 

V  i-«s  /=  #¥. 

7 — ?  /yv/  \  > 

7  ’  V>/  /-  -«-■* 


^  /  *  t  1  *  K  ^  -o 

/  ; !  -K 

-  {,;•;  M  vs 

\r :  j  ~f  \  / 

^  /  *  /  /-  3  5 

/  y  1  /'  ? 

v 

xw:m?rj:-- 


•7?ri?  Vl 


v~  ~  vvl' 

MS 


Srr 


c.  - 


<D 

*  Jh  4 

^ !  j* 

i  I  § 

r  «  a) 

H  Q 


■;  "M-/ 

\w  >/ 


*A  * 

■  M  M":--  ? 

Z  -\  O  * 

-  1/  8 
*>  \  s  V-* 


S 


:  /^i 


^7, 

<<  5 


Z  ^v. /fN  , 


vm 


-  ^ 

iM. 


=\i  % 

\  »  _JnT 


X  '5 


8  ,18 


^7  ^  ^ 

s '  /a  >~2 


Fig.  2-7.  Ballistic  temperature  and  density  analysis  using  CRAM:  0600  MST,  25  January  1965  (percent 
departure  from  Standard). 
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Fig.  2-11.  Ballistic  temperature  and  density  analysis  using  CRAM:  1400  MST,  25  January  1965  (percent 
departure  from  Standard). 


points.  Other  degrees  of  smoothing  will  be  tried  in  the  future. 

Analyses  wei  ^produced  for  each  of  the  ten  artillery  zones  using  the  procedures 
outlined  above.  These  zone  analyses  were  then  combined  into  ballistic  analyses  by 
computing  the  weighted  sums  using  the  weights  given  in  Table  2-i.  Figures  2-2 
through  2-6  show  the  five  analyzedfields  of  ballistic  u-  and  v-wind  components  (in 
knots)  for  0600,  0800,  1000,  1200,  and  1400  MST.  The  temperaiure  and  density  analyses 
are  given  in  Figs.  2-7  through  2-11  (departure  from  Standard,  in  percent).  The  bal¬ 
listic  wind  analyses  appear  much  different  than  those  produced  by  hand  in  the  displaced 
stations  experiment  (Figs.  2-7  through  2-11  in  Quarterly  Report  No.  2).  This  is  no 
doubt  due  to  the  fact  that  no  initial  guess  is  used  in  the  hand  analyses  which  leads  to  a 
considerable  amount  of  subjectivity  in  areas  of  map  where  there  are  no  data. 

The  u-component  ballistic  wind  analyses  show  a  gradual  weakening  trend 
through  the  entire  map  series.  A  minimum  of  84  knots  northwest  of  Douglas  can  be 
followed  westward  map  by  map  reaching  an  area  just  west  of  Ft.  Huachuca  at  1400  MST. 
More  testing  may  be  required  to  see  if  this  pattern  is  realistic  or  whether  it  is  due  to 
some  characteristic  of  the  analysis  technique. 

Small  negative  values  of  the  v-component  of  ballistic  wind  are  present  over  the 
southeastern  part  of  the  analysis  area  at  0600  and  0800  MST.  At  1000  MST  there  is  a 
rather  sudden  change  in  the  eastern  part  of  the  region  to  large  negative  values  which 
remain  for  the  duration  of  the  series. 

The  temperature  and  density  analyses  (Figs.  2-7  through  2-11)  show  relatively 
small  gradients  with  the  range  of  values  over  an  entire  map  being  generally  less  than 
1%  of  Standard.  The  temperature  values  show  a  cooling  trend  through  the  series, 
reflecting  a  cooling  trend  in  the  upper  zones.  The  surface  heating  affecting  the  lower 
zones  does  not  contribute  much  to  the  overall  ballistic  temperatui;e  because  of  the  way 
the  individual  zones  are  weighted  (see  Table  2-1). 

The  density  fields  exhibit  increasing  gradients  during  the  map  series.  The  maxi¬ 
mum  values  are  located  in  the  southern  part  of  the  analysis  area,  with  a  ridge  extend¬ 
ing  northeastward  and  becoming  more  pronounced  with  time. 

The  ballistic  quantities  of  u,  v,  temperature,  and  density  were  computed  from 
these  analyses  at  the  approximate  midpoint  of  the  artillery  trajectories  and  are  shown  in 
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Table  2-2.  Also  shown  are  those  quantities  computed  from  the  hand-produced  map 
series  contained  in  Quarterly  Report  No.  2.  In  using  CRAM,  the  changes  of  the  values 
with  time  in  the  firing  area  are  generally  more  gradual  than  the  changes  shown  by  the 
other  analysis  methods  used.  This  is  due,  at  least  in  part,  to  the  incorporation  of  time 
in  the  generation  of  the  initial-guess  fields. 

Of  particular  interest  is  the  application  of  these  ballistic  corrections  to  the 
actual  artillery  firings  that  took  place  at  1000  and  1400  MST  for  the  same  day.  The 
provided  gun  data  and  unit  effects  were  employed  for  this  purpose.  Ballistic  tempera¬ 
ture  and  density’  values  were  needed  to  compute  the  overall  correction.  These  were 
computed  from  the  CRAM  analyses  (Figs.  2-9  and  2-11}  with  the  values  of  0.34  for 
temperature  and  0.13  for  density  at  1000  MST  and  -0.02  for  temperature  and  0.25  for 
density  at  1400  MST.  From  fixed  station  analyses  (FXD-1),  the  values  were  calculated 
to  be  0.15  and  0.64  at  1000  MST  and  -0.08  and  0.75  at  1400  MST.  Because  no  values  of 
temperature  and  density  were  available  from  the  DISPL  and  FXD-2  experiments,  the 
FXD-1  values  were  used  for  FXD-2  and  the  CRAM  values  were  used  for  DISPL.  The  results 
are  shown  in  Table  2-3.  The  target  (hypothetical)  coordinates  are  given  in  meters  in 


TABLE  2-2 

COMPARISON  OF  BALLISTIC  WIND  COMPONENTS  AT 
MIDPOINT  OF  PROJECTILE  TRAJECTORY* 


Time 

u- 

-component  (knots] 

v- 

-component  (knots) 

(MST) 

FXD-1 

DISPL 

FXD-2 

CRAM 

FXD-1 

DISPL 

FXD-2 

CRAM 

0600 

86.2 

89.7 

91.0 

90.8 

— 

-  8.5 

-10.7 

0800 

82.0 

77.5 

84.0 

81.3 

-  9.2 

-16.5 

-12.0 

1000 

72.8 

69.2 

*5.0 

74.8 

-16.5 

-16.7 

-21.5 

-16.0 

1200 

70.7 

60.3 

75.0 

68.1 

-17.4 

-22.2 

-17.8 

1400 

64.8 

52.2 

68.0 

61.0 

-17.6 

-17.5 

-15.8 

*  FXD-1:  hand  analysis  based  on  fixed  stations. 

FXD-2:  hand  analysis  based  on  fixed  stations,  re-analyzed. 
DISPL:  hand  analysis  based  on  radiosonde  displacement. 
CRAM:  analysis  based  on  CRAM  objective  analysis. 
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the  Universal  Transverse  Mercator  grid  system.  The  range,  deflection,  and  vector 
errors  designated  “NO  CORRECTIONS"  are  based  on  the  distance  between  the  target 
and  the  center  of  impact.  The  range,  deflection,  and  vector  errors  -WITH  CORREC¬ 
TIONS"  refer  to  the  corrected  impact  point  after  the  drift,  coriolis,  and  meteorological 
effects  have  been  accounted  for.  (Negative  range  errors  refer  to  impact  short  of  the 
target  and  negative  deflection  errors  refer  to  impact  to  the  left  of  the  target  when 
viewed  from  the  gun  location.)  Also  given  are  the  x  and  y  components  of  the  vector 
error  which,  when  added  to  the  target  coordinates,  give  the  coordinates  of  the  cor¬ 
rected  impact  point.  The  corrections  due  to  temperature,  densit}*,  and  range  wind 
effect  only  the  range  (negative  corrections  decrease  the  range),  while  corrections  due 
to  cross  wind  effect  only  the  deflection  (negative  corrections  displace  the  impact  point 
to  the  left).  The  variability  of  the  residual  errors  among  the  several  approaches  tested 
is  due  inainty  to  the  differences  in  ballistic  winds,  with  the  differences  due  to  temper¬ 
ature  and  density  playing  a  relatively  minor  role.  These  residuals  are  plotted  relative 
to  the  target  in  Fig.  2-12.  In  most  of  these  impacts,  the  range  error  is  negative  for  the 
corrected  impact  points.  It  will  be  shown  in  the  next  section  that  this  negative  range 
error  was  rather  2Jiamoious  when  compared  to  several  other  firing  situations. 
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TEST  OF  METEOROLOGICAL  EFFECTS 
(25  January  1965,  Gun  No.  1,  AM  &  PM) 
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TABLE  2-3  (continued) 
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Fig.  2-12.  Residual  errors  for  1000  MST  (AM)  and  1400  MST  (EM)  firings  (CM  No.  1)  on  25  January  1965 
using  corrections  based  on  various  analysis  procedures  (from  Table  2-3). 


3.0  CORRECTIONS  BASED  ON  A  SINGLE  STATION 


A  small  computer  program  was  written  for  the  IBM  1620  which  makes  the 
necessary  geometrical  and  trigonometrical  steps  to  systematically  modify  a  given 
impact  point  with  given  meteorological  parameters  and  unit  effects  in  order  to  yield  a 
corrected  impact  point.  The  program  prints  out  the  corrections  and  the  various  error 
statistics.  By  computing  this  information  for  all  available  firing  times  and  using  the 
Ft.  Huachuca  sounding  as  the  only  basis  for  computing  meteorological  questions,  the 
resulting  data  could  be  used  to  study  the  general  pattern  of  the  impact  points.  The 
data  also  provided  a  useful  control  in  comparing  results  from  the  objective  analysis 
technique. 

The  gun  data  resulted  from  the  firings  of  two  8-inch  howitzers,  each  firing  about 
ten  rounds,  alternatively,  at  one  minute  intervals.  For  a  given  day  there  were  two 
firing  times:  1000  MST  and  1400  MST.  The  locations  of  the  individual  impact  points 
were  averaged  to  give  a  center  of  impact  (Cl)  for  each  gun’s  firing.  These  data  were 
supplied  by  the  sponsor.  The  coordinates  of  the  two  guns  and  their  azimuths  are  as 
follows: 

Gun  No.  1:  3493987.7N  555472.8E  1325.2  mils,  Grid  North 
Gun  No.  2:  3494002.5N  555469.4E  1322.8  mils,  Grid  North 

Of  the  24  firing  times,  there  were  19  for  which  the  Ft.  Huachuca  sounding 
yielded  line-10  ballistic  data.  The  corrections  were  applied  to  the  appropriate  Cl’s  by 
using  the  IBM  1620  program  written  for  that  purpose.  The  results  are  given  in 
Table  3-1  and  in  Figs.  3-1  and  3-2,  which  show  the  corrected  Cl’s  relative  to  the  tar¬ 
get.  It  is  immediately  obvious  that  there  is  a  preponderance  of  corrected  impact  points 
for  which  the  range  error  is  positive  (too  long).  The  reason  for  this  is  not  known  at 
the  present  t  me.  There  does  not  appear  to  be  any  noticeable  bias  in  the  deflection 
errors.  Table  3-2  shows  a  summary  of  the  error  statistics.  In  addition  to  the  range 
bias,  the  total  shows  that  there  is  no  appreciable  difference  in  errors  resulting  from 
Gun  No.  1  and  Gun  No.  2  firings. 
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*Ft.  Huachuca  soundings  were  not  available  for  Cases  11, 12, 13, 14, 17, 18, 31, 32, 41,  and  42 
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TABLE  3-1  (continued) 
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Residual  errors  for  Gun  No.  1  firings  listed  in  Table  3-1.  Numbers  refer  to  case  numbers  in 


Fig.  3-2.  Residual  errors  for  Gun  No.  2  firings  listed  in  Table  3-1.  Numbers  refer  to  case  numbers  in 


TABLE  3-2 

SUMMARY  OF  RESIDUAL  ERRORS  FROM  TABLE  3-1  FOR 
19  FIRING  TIMES  (in  m) 


Gun 

no. 

Average  error 
(algebraic) 

Average  error 
(absolute) 

Residual 

vector 

error 

Location  of 
impacts 

Range 

Deflection 

Range 

Deflection 

Long 

Short 

.  . 

Left 

Right 

1 

-5-187 

+  3 

244 

85 

272 

14 

5 

8 

11 

2 

+198 

-29 

263 

88 

287 

14 

*  5 

10 

9 

4.0  POST -PROCESSING  PROGRAM 

A  computer  program  is  being  written  for  the  IBM  7090  which  will  be  an  expanded 
version  of  the  small  IBM  1620  program  mentioned  in  the  previous  section.  It  will  com¬ 
pute  fields  of  residual  errors  which  can  be  averaged  over  a  series  of  maps  to  deter¬ 
mine  spatial  and  temporal  variability  of  artillery  corrections.  Corrections  based  on 
concurrent  meteorological  effects  will  be  compared  with  corrections  based  on  such 
effects  2-,  4-,  6-,  and  8-hours  old.  The  change  in  ballistic  quantities  from  0600  to 
0800  MST  will  be  used  to  generate  a  predicted  1000  MST  set  of  ballistic  maps  which 
can  be  used  to  compute  corrections.  Likewise,  a  1400  MST  prediction  can  be  gener¬ 
ated  from  the  change  between  0600  and  1000  MST  as  well  as  the  change  between  1000 
and  1200  MST. 
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